2009 Fermi Symposium, Washington, B.C., Nov. 2-5 



1 



Extending the Galactic Cosmic Ray electron + positron spectrum measured by 

the Fermi LAT 
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Launched on the 11th of June 2008, the Fermi Large Area Telescope (LAT) has made several 
outstanding scientific contributions to the high energy astrophysics community. One of these contri- 
butions was the high statistics measurement of the Galactic Cosmic Ray (GCR) electron + positron 
spectrum from 20 GeV to 1 TeV. The Fermi satellite is in a nearly circular orbit with an inclination 
of 25.6 degrees at an altitude of 565 km. Given this orbit it is possible to measure the GCR electrons 
+ positrons down to roughly 5 GeV. However, this lower limit in energy is highly dependent on 
the orbital position of the LAT in geomagnetic coordinates due to the rigidity cutoff. In order to 
measure the spectrum down to these energies it is necessary to sample the population of electrons + 
positrons in several different geomagnetic positions. In this poster we present the analysis performed 
to extend the lower limit in energy of the GCR electron + positron spectrum measured by the Fermi 
LAT. 



I. INTRODUCTION 

The electron component of the Galactic Cosmic 
Ray (GCR) radiation is widely recognized as a unique 
probe to address a number of significant questions 
concerning the origin of cosmic rays and their propa- 
gation in our galaxy (see @ for a review). For energies 
above a few GeV the dominant energy loss mecha- 
nisms for electrons and positrons are the synchrotron 
and inverse compton processes. These processes have 
an energy loss rate whose magnitude increases with 
the square of the electron energy and as a consequence 
these electrons have a radiative lifetime inversely pro- 
portional to their energy Q. Because of this rela- 
tively short radiative lifetime, the GCR electrons can- 
not travel intergalactic distances through the cosmic 
microwave background and thus excludes an extra- 
galatic contribution to the measured spectra. After 
the first six months of nominal science data taking, 
the Fermi LAT was able to make the first systematics 
limited measurement of the GCR e~ + e + (referred to 
as electrons for the rest of this work) spectrum from 
20 GeV to 1 TeV 2]. This measurement together 
with the recently published results from experiments 
such as PAMELA 0], ATIC Q and HESS @ has set 
the stage for new and exciting GCR science. To help 
make constraints on possible propagation models it 
is important to extend the measurement of the GCR 
electrons to lower energies. This task is the main topic 
of this work and the analysis used to extend the GCR 
electron spectrum measured by the Fermi LAT will be 
described in detail in this paper. 



II. INSTRUMENT DESCRIPTION 

The Large Area Telescope, the main instrument on- 
board the Fermi observatory, is a pair-conversion tele- 
scope composed of: 

(i) a tracker-converter, including 18 Silicon Strip 



Detector (SSD) x—y tracking planes, interleaved 
with tungsten conversion foils; 

(ii) a Csl imaging calorimeter; 

(iii) a segmented Anti- Coincidence Detector (ACD) 
surrounding the tracker subsystem. 

Though it essentially follows in the footsteps of the 
previously flown pair-conversion telescopes (particu- 
larly EGRET P|) in its basic design, the LAT is 
equipped with the state of the art in high-energy 
physics detector technology and features many design 
improvements with respect to its predecessor, result- 
ing in a significant leap in sensitivity (for a detailed 
description of the detector see @). 

From the standpoint of the electron detection, the 
tracker section effectively acts as a finely segmented 
pre-shower (1.5 radiation lengths on axis) which allows 
to image the initial part of the electromagnetic show- 
ers, providing a crucial handle for the background re- 
jection. The hodoscopic configuration of the calorime- 
ter (8.6 radiation lengths on axis) provides a full 3- 
dimensional reconstruction of the shower shape which 
is essential both for the energy reconstruction and 
for the discrimination between electromagnetic and 
hadronic showers. The ACD also provides a relevant 
contribution to the background rejection. 



A. Trigger and filter 

Two of the main design improvement of the LAT 
are a flexible trigger logic and onboard event filter- 
ing. The chosen detector technology allows a mini- 
mum dead time per event of about 26 /is (essentially 
due to instrument readout). Therefore the LAT can 
trigger at the passage of (almost) every particle, in- 
cluding cosmic rays (only events that clearly behave 
as MIPs or heavy ions are hardware-prescaled to avoid 
the induction of unnecessary deadtime). This is a 
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revolution with respect to the previous generation of 
pair conversion telescopes, in which the instrumental 
deadtime forced to apply a harsh photon selection at 
the level of the hardware trigger. The LAT trigger 
rate is of the order of few kHz but due to telemetry 
bandwidth limitations an onboard event selection is 
applied and roughly half of the events are actually 
written on spacecraft memory. This is accomplished 
by several software configurable onboard filters. The 
rate of events downlinked to Earth is of the order of 
few hundred Hz, still mostly charged particles (i. e. 
protons, electrons, GCR metals). The refined anal- 
ysis, both for the photon (cfr.@) and CR electron 
science is performed on the ground. For more details 
about the on-orbit rates, see Q. In the routine data 
taking configuration the onboard filter instance pro- 
viding the main data source for photon science (the 
GAMMA filter) is configured to pass all the events de- 
positing more than 20 GeV in the calorimeter. There- 
fore it constitutes the main source of data for the high- 
energy electron analysis, as well. In addition to that, 
an unbiased sample of all trigger types, prescaled by 
a factor of 250, is also downlinked for diagnostic pur- 
poses (through the so called DGN filter) ; it is perfectly 
adequate to study the electrons at relatively low ener- 
gies, where the higher fluxes compensate for the lower 
statistics due to the event prescaling. The events col- 
lected through the DGN filter are essential for the ex- 
tension of the GCR electron spectrum, in fact they 
constitute the main source of electrons for this study. 



III. GEOMAGNETIC ENVIRONMENT 

The Earth's field is closely approximated by the 
field of a dipole positioned near the center of the 
Earth. For a dipole field, 

R = Rq cos 2 (1) 

defines lines of constant magnetic field, also known as 
"L shells". Ro is the radial distance to the field line 
where it crosses the geomagnetic equator and R is the 
radial distance to the point where the field is B at lat- 
itude 9. From equation Q] it is possible to define the 
coordinate L, which is defined as L = Ro/Re where 
Re is the Earths's radius (6371 km) [12J. Therefore, 
positions around the Earth with the same value of 
L are magnetically equivalent. Due to the fact that 
the Earth's dipole is offset and tilted with respect to 
the center of the Earth, L values must be calculated 
based on a detailed field model [l[ and are a function 
of geomagnetic latitude [Hj]. The shielding effect of 
the Earth's magnetic field is smaller (larger) for larger 
(smaller) values of L and with Fermi's nearly circu- 
lar orbit it is possible to measure the GCR electron 
spectrum in the range 1.0 < L < 1.73 which trans- 
lates into a minimum rigidity of ~ 6 GeV. Figure [T] 



illustrates the distribution of L and the correspond- 
ing vertical rigidities for the Fermi orbit, based on the 
IGRF model [l(. It important to stress here that the 
vertical cutoff rigidity distribution shown in figure [T] 
is simply an illustration and we do not depend on this 
model in any way for this analysis. 

Fermi's orbit intersects the South Atlantic Anomaly 
(SAA) which is a region of the Earth's radiation belt 
which features geomagnetically trapped protons with 
energies up to hundreds of MeV and electrons with 
energies up to tens of MeV. The extreme conditions 
within this region impose constraints on the LAT op- 
erations, in particular the triggering, recording and 
transmission of science data are stopped during the 
SAA passages Q. The SAA is represented by the 
dark shaded region in figure [TJ 

IV. EVENT SELECTION 

The LAT is sensitive to electrons over more than 
five orders of magnitude in energy (~ 100 MeV - ~ 1 
TeV). Across this huge energy range both the typical 
event topology and the relative fractions of signal and 
background in the cosmic-ray flux impinging on the 
detector undergo dramatic variations. For these rea- 
sons it's not trivial to develop a single, unitary analy- 
sis strategy providing the necessary electron detection 
efficiency and the hadron rejection power across the 
whole energy range. In fact we have elaborated two in- 
dependent electron selections, tuned for relatively low 
energies and high energies respectively (which we shall 
refer to as LE and HE in the rest of the paper). The 
split point, in energy, between these two analyses can 
be naturally placed at 20 GeV for at least two reasons. 
The first one is that the onboard filtering is essentially 
disengaged above this energy (whereas we only have a 
prescaled unbiased sample below 20 GeV). The other 
one is that below k, 20 GeV, the shielding effect of the 
geomagnetic field (which depends on the geomagnetic 
latitude and, hence, on the position of the spacecraft 
across its orbit) becomes important. In fact the he se- 
lection has been mostly developed exploiting the data 
source provided by the GAMMA filter for the purpose 
of measuring the GCR electrons, whereas the le se- 
lection was designed to measure both the GCR and 
the electron albedo population, using the DGN filter 
data source. 

It must be stressed, however, that these divisions 
are partially artificial and there's a significant overlap 
in energy between the two selections. This can be used 
as an independent cross check for both of the analyses. 
The details of the selection criteria for the le selection 
can be found in [l(| and for the HE in [2j. 

One of the important figures of merit which come 
out of the selection is the effective geometry factor, 
defined as the effective area integrated over the field 
of view of the detector, and is shown in figure O The 



2009 Fermi Symposium, Washington, B.C., Nov. 2-5 



3 



McllwainL 




FIG. 1: Map of Mcllwain L values for the Fermi orbit. Overlaid in contours are the corresponding values for vertical 
rigidity cutoff. These values were calculated using the 10 th generation IGRF model. It is important to note that the 
model is not valid inside the South Atlantic Anomaly (SAA) region, illustrated by the dark gray polygon in this figure. 
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FIG. 2: Effective Geometric Factor for LE (blue squares, 
multiplied by a factor 250) and he events (red triangles). 
The results for the LE selection are still preliminary. 



LE geometry factor has multiplied by the DGN filter 
prescale factor of 250 for graphical clarity and in or- 
der to give an idea of the relative electron efficien- 
cies of the two selections. The shape of the geome- 
try factor is almost entirely due to the event selec- 
tion that must obtain a good background rejection 
power in the whole energy range. The sharp decrease 
in HE selection just below 30 GeV is mainly due to 
the GAMMA filter onboard event filtering designed to 
remove charged particles with deposited energy lower 
that 20 GeV. 



V. SPECTRUM RECONSTRUCTION 

The electron flux is evaluated by subtracting the 
residual background from the count rate and then cor- 
recting the observed number of events with the geom- 
etry factor. The evaluation of the residual hadron 
background takes advantage of the detailed on-orbit 
particle environment model built by the LAT collab- 



oration. This highly detailed model has been inten- 
sively used to develop all the 7-ray background re- 
jection algorithms, both on-orbit and off-line. The 
model includes cosmic rays and earth albedo 7-rays 
starting from 10 MeV and is valid outside the South 
Atlantic Anomaly (since the LAT does not take data 
while passing through the SAA). The particle fluxes 
are chosen to fit experimental data of several past ex- 
periments, more details can be found in [||. 

The event selection is tuned to keep the residual 
contamination below about 20% in the whole energy 
range. The sample of events passing the GAMMA filter 
are the main source of GCR electrons given that these 
are all those events with raw energy > 20 GeV and 
for these energies the effects of the Earth's magnetic 
field on the incoming GCR electron flux are negligible. 
In order to extend the energy range of the measured 
spectrum (down to lower energies) it is necessary to 
use the events passing the DGN filter. This extension 
is not a straight forward task and requires to take ad- 
vantage of the fact that the rigidity cutoff varies as a 
function of orbital position and consequently measure 
the spectra in various bins of geomagnetic position. 
As already mentioned in section [TJ the shielding ef- 
fect of the Earth's magnetic field is smaller for larger 
values of the L coordinate and therefore to measure 
electrons of Galactic origin with energies in the lower 
end of our rigidity cutoff range it is necessary to sam- 
ple all those events collected when the detector was 
located at large values of L. Each L bin has an associ- 
ated rigidity cutoff value and it is possible to measure 
this value by considering that the shape of the primary 
spectrum around the geomagnetic rigidity cutoff can 
be parametrized as: 



ME) 



1 



(2) 



1 + (E/E C )-* 
Therefore the full spectrum can be fitted with a func- 
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FIG. 3: Left panel: The measured electron flux in three Mcllwain L bins. For each bin the fit of the flux with equation 
[3] and the resulting estimated cutoff rigidity, E c , is shown. As described in the text, E c decreases for larger values of 
Mcllwain L. Right panel: Illustration of the reconstruction technique used to measure the the primary cosmic-ray electron 
spectrum down the lowest accessible energy for the Fermi orbit. The blue points represent the GCRE flux multiplied by 
the energy cubed obtained via the le analysis and the red points are from the he analysis. The solid lines depict the 
systematic uncertainties and the vertical dashed lines serve to show from which Mcllwain L region of the orbit the flux 
was measured. These results are preliminary. 



tion of the form: 

dE- CsE + 1 + (E/E C )S (3) 

Where the subscript s stands for secondary and repre- 
sents the albedo population of electrons while the sub- 
script p stands for primary and represents the Galactic 
component of the spectra and T is the spectral index. 
The value of E c , obtained by the fit is an estimation of 
the rigidity cutoff. An example of such a fit is shown 
in figure [3] for three L bins. Due to the complexity 
of particle orbits in the Earth's magnetosphere, the 
geomagnetic cut-off is not sharp but rather a smooth 
transition and therefore it is necessary to include a 
multiplicative factor, X, to the value of E c obtained 
by the fit. This factor serves as a sort of padding, and 
its value is driven by the need to balance between two 
requirements, namely to measure the lowest possible 
Galactic CRE energy while at the same time making 
sure that the measured flux value is not affected by the 
shielding due to the Earth's magnetic field. The value 
of X which satisfies these requirements was chosen by 
performing a scan on several values of X around unity. 
Having done this we are able to quantify the changes 
in the flux level induced by the value of X used to 
extend the spectrum. From this study we found that 
a value of 1.15 satisfies the above requirements. 

Once the value for E c has been estimated for each 
L bin and the padding size has been fixed, it is pos- 
sible to extend the GCRE spectrum. An example of 
how this is done is shown in the right panel of figure 
[3] where the vertical dashed lines denote the L slice 
from which the flux was measured. Where L>1.00 
corresponds to the orbital averaged data (i.e. where 
the shielding effects from the Earth's magnetic field 
are negligible) and L>1.70 is the highest value for L 
sampled by the Fermi orbit. The spectra obtained 



with le and the he are shown in this figure to illus- 
trate the very good agreement in the overlap region. 
The solid (blue and red) lines represent the systematic 
uncertainties for the two selections. 



VI. DISCUSSION 

The GCR electron spectrum measured by the Fermi 
LAT from - 7 GeV to 1 TeV is shown in figured to- 
gether with the measurements made from ground, bal- 
loon and spaced based experiments. The overlap re- 
gion between the le and the he selections for energies 
between 20 and 80 GeV shows very good agreement 
and serves as an excellent cross check between the two 
separate analyses. The low energy part of the spec- 
trum measured by Fermi shows the same rising trend 
in flux as both the AMS and HEAT measurements 
and is compatible with those measurements within the 
systematic errors. The extension of the GCR elec- 
tron spectrum provides valuable information needed 
to make constraints on possible propagation models 
of the GCRs. However, a discussion on the interpre- 
tation of this new measurement is beyond the scope 
of this paper. 
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FIG. 4: The GCR electron spectrum measured by the Fermi LAT, red points correspond to the he selection and blue 
points to the le. For comparison the measurements from the experiments AMS, HEAT, ATIC and HESS are also shown. 
The extension of the Fermi LAT measurement is in good agreement within the systematic errors with the measurements 
made by AMS and HEAT. 
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